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ABSTRACT 

In several merging clusters of galaxies so-called cluster radio relics have been observed. 
These are extended radio sources which do not seem to be associated with any ra- 
dio galaxies. Two competing physical mechanisms to accelerate the radio emitting 
electrons have been proposed: (i) diffusive shock acceleration and (ii) adiabatic com- 
pression of fossil radio plasma by merger shock waves. Here the second scenario is 
investigated. We present detailed 3-dimensional magneto-hydrodynamical simulations 
of the passage of a radio plasma cocoon filled with turbulent magnetic fields through a 
shock wave. Taking into account synchrotron, inverse Compton and adiabatic energy 
losses and gains we evolved the relativistic electron population to produce synthetic 
polarisation radio maps. On contact with the shock wave the radio cocoons are first 
compressed and finally torn into filamentary structures, as is observed in several clus- 
ter radio relics. In the synthetic radio maps the electric polarisation vectors are mostly 
perpendicular to the filamentary radio structures. If the magnetic field inside the co- 
coon is not too strong, the initially spherical radio cocoon is transformed into a torus 
after the passage of the shock wave. Very recent, high-resolution radio maps of cluster 
radio relics seem to exhibit such toroidal geometries in some cases. This supports the 
hypothesis that cluster radio relics are fossil radio cocoons that have been revived by 
a shock wave. For a late-stage relic the ratio of its global diameter to the filament di- 
ameter should correlate with the shock strength. Finally, we argue that the total radio 
polarisation of radio relic should be well correlated with the 3-dimcnsional orientation 
of the shock wave that produced the relic. 

Key words: galaxies: clusters: general - intergalactic medium - shock waves - mag- 
netohydrodynamics - radio continuum: general - polarization 



1 INTRODUCTION 

In the current picture of hierarchical structure formation 
clusters of galaxies grow mainly by the merging of smaller 
and moderately sized sub-clusters. During such merger 
events a significant fraction of the kinetic energy of the in- 
ter galactic medium (IGM) is dissipated in Mpc-sized shock 
waves. The shock waves are responsible for heating the IGM 
to temperatures of several keV. Moreover, they are likely 
to inject and accelerate relativistic particle populations on 
cluster scales. 

Cluster- wide relativistic electron populations are indeed 
observed in several merging or post-merging clusters. These 
electrons reveal their presence by emitting synchrotron emis- 
sion in the magnetic fields of the IGM to form the so-called 
cluster radio halos. Cluster radio halos have steep radio spec- 
tra (spectral index a m 1 — 1.5); they are unpolarised, and 
their diffuse morphologies are often similar to those of the 
thermal X-ray emission of the cluster gas (Govoni et al. 
2001). A comparison of the typical radiative lifetime of the 



radio emitting electrons, which is of the order of 0.1 Gyr, to 
the typical dynamical timescale of the cluster merger, which 
is of the order of 1 Gyr, shows that the electron population 
of a halo cannot be produced directly by a merger shock. Ei- 
ther these electrons have been re-accelerated more recently, 
or they have been injected by a long-lived shock-accelerated 
proton population via the hadronic charged pion produc- 
tion. For reviews on recent observations of diffuse cluster 
radio sources see Feretti (1999), Giovannini, Tordi, Feretti 
(1999), and Kempner & Sarazin (2001) . 

Unlike the radio halos, the so-called cluster radio relief 



The name 'cluster radio relic ' was originally chosen in order to 
reflect the similarities of these sources with 'relic radio galaxies'. 
The latter are dying radio cocoons that are no longer supplied 
with fresh radio plasma from the central parent galaxy. In clus- 
ter relics, usually no parent galaxy can be identified. Moreover, 
they have only been found in clusters of galaxies, and therefore 
they were named 'cluster radio relics'. Recently, the connection 
between cluster radio relics and radio galaxies has been debated. 
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Figure 1. Cluster radio relic in Abell 85 at 327 MHz (grey-scale, 
Giovannini & Feretti 2000) and at 1.4 GHz (contours, Slee et al. 
2001). The 327 MHz maps demonstrates that the thin 1.4 GHz 
filament extending Northwards belongs in fact to a closed loop in 
projection. The two smaller, southern blobs on the 1.4 GHz map 
have also torus-like morphologies. 

are more directly related to cluster mergers (EnBlin et al. 
1998, Roettiger et al. 1999, Venturi et al. 1999, Enfilin & 
Gopal-Krishna 2001). They are also extended radio sources 
with a steep spectrum. In the literature radio relics are of- 
ten confused with radio halos even though several distinctive 
properties exist. Cluster radio relics are typically located 
near the periphery of the cluster; they often exhibit sharp 
emission edges and many of them show strong radio polari- 
sation. 

In several cases it could be shown that shock waves are 
present at the locations of the relics. In Abell 2256 and Abell 
1367 temperature substructures of the hot IGM could be 
detected (Briel & Henry 1994, Donnelly et al. 1998) which 
points towards the presence of shock waves at the location of 
the cluster relics in these clusters. For Abell 754 (Roettiger 
at al. 1998; Kassim et al. 2001), Abell 2256 (Roettiger et 
al. 1995), Abell 3667 (Roettiger et al. 1999) and also the 
Coma cluster (Burns et al. 1994) numerical simulations of 
mergers were fitted to the X-ray data. These simulations 
predict shock waves at locations of observed cluster radio 
relics. 

The cluster radio relic 1253+275 in the Coma cluster 
shows a morphological connection to the nearby radio galaxy 
NGC 4789 (Giovannini, Feretti & Stanghellini 1991). This 
suggests that radio relics may be fossil radio plasma that 
has been revived by a shock. Fossil radio plasma is the for- 
mer outflow of a radio galaxy in which the high-energy radio 

Alternative, possibly less misleading names have been proposed, 
e.g. 'radio flotsam' by R. Ekers. Here, we use the original termi- 
nology a) in order to be consistent with the former literature and 
b) since our work provides support for a scenario, in which these 
objects consist of re-illuminated fossil radio plasma, so that they 
are indeed shock processed relic radio galaxies. 



emitting electrons have lost their energy. Due to their invis- 
ibility in the radio these cocoons are also called radio ghosts 
(Enfilin 1999). 

The first relic formation models considered diffusive 
shock acceleration (Fermi I) as the process producing the 
radio emitting electrons (Enfilin et al. 1998, Roettiger et 
al. 1999, Venturi et al. 1999). However, when a fossil ra- 
dio cocoon is passed by a cluster merger shock wave, with 
a typical velocity of a few 1000 km/s, the cocoon is com- 
pressed adiabatically and not shocked because of the much 
higher sound speed within it. Therefore, shock acceleration 
cannot be the mechanism that re-energises the relativistic 
electron population. But the energy gained during the adi- 
abatic compression combined with the increase in the mag- 
netic fields strength can cause the fossil radio cocoon to 
emit radio waves again. One prerequisite for this is that the 
electron population is not older than 0.2-2 Gyr (Enfilin 
& Gopal-Krishna 2001). The timescale depends on the con- 
ditions in the surroundings, mainly the external pressure. 
In high-pressure environments, such as in cluster cores, the 
synchrotron losses are expected to be much higher due to 
the higher internal magnetic fields of pressure-confined ra- 
dio plasma. This leads to a shorter maximum age for the 
fossil radio plasma if it is to be revived. Therefore, relics 
are found preferentially at the periphery of clusters where 
the pressure is lower. Moreover, numerical simulations show 
that merger shocks are found more frequently in peripheral 
cluster regions than in the denser cores (Quilis, Ibanez, & 
Saez 1998; Miniati et al. 2000). Both effects, the longer time 
for radio plasma to be revivable, and the higher frequency of 
shock waves could explain why cluster radio relics are more 
frequently observed at the outskirts of clusters than in more 
central regions. 

Enfilin & Gopal-Krishna (2001) showed that the spec- 
tral properties of cluster radio relics are well reproduced 
by this scenario. Here, we demonstrate that the observed 
morphologies and polarisation patterns are reproduced by 
this model as well. This is done with the help of the first 
3-dimensional magneto-hydrodynamical (MHD) simulations 
of a fossil radio cocoon that is passed by a shock wave. We 
produce artifical radio maps that can be compared directly 
to high-resolution radio maps of relics. 



2 METHOD 

2.1 3D MHD Simulation 

The magneto-hydrodynamical simulations were obtained us- 
ing the ZEUS-3D code which was developed especially for 
problems in astrophysical hydrodynamics (Stone & Norman 
1992a, b). The code uses finite differencing on a Eulerian grid 
and is fully explicit in time. It is based on an operator-split 
scheme with piecewise linear functions for the fundamental 
variables. The fluid is advected through a mesh using the up- 
wind, monotonic interpolation scheme of van Leer. The mag- 
netic field is evolved using a modified constrained transport 
technique which ensures that the field remains divergence- 
free to machine precision. The electromotive forces are com- 
puted via upwind differencing along Alfven characteristics. 
For a detailed description of the algorithms and their nu- 
merical implementation see Stone & Norman (1992a, b). 
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Figure 2. Evolution (from the initial configuration (left) to the final one (right)) of the gas density (top) and magnetic field energy 
density (bottom) in a central slice through the simulation volume. The pre-shock region is on the left hand side, and the post-shock 
region on the right hand side of the stationary shock wave, which is located at the centre of the computational box. The formation of a 
torus is best visible in the density plot. In this simulation the projection of the vector potential of the magnetic fields is not done by Eq. 

but by the modified form A = [I — ff(r)(I — ff)]/(l — 2g(r)/3) ■ A in order to have a more uniform magnetic field within the radio 
cocoon. The numbers on the axes label the zones of the grid. 



The simulations were computed on a Cartesian grid 
with 100 3 and 200 3 equally spaced zones. The fluid flows 
in the :r-direction with inflow boundary conditions at the 
lower boundary and outflow conditions at the outer bound- 
ary. The boundary conditions in the y- and ^-directions were 
chosen to be reflecting. The simulation was set up such that 
a stationary shock formed at the centre of the simulation 
box that is perpendicular to the direction of the flow. 

In the pre-shock region a spherical bubble of radius R 
was set up, in which the density was lowered by a factor of 10 
with respect to the environment. In turn, the temperature 
in the bubble was raised such that the bubble remained in 
pressure equilibrium with its surroundings. Inside the bub- 
ble a magnetic field was set up as follows: The magnetic 
vector potential A was composed of 125 Fourier modes with 
random phases and amplitudes. A was then scaled as fc~ 3 , 
where k is the wavenumber of the mode. This leads to a 
magnetic power spectrum of Ps(fc)dfc ~ k~ l dk. To avoid 
leakage of the magnetic field into the surroundings of the 
bubble, we employed the projection: 

A = [I- ff (r)(I-ff)]-A, (1) 

where I is the identity, r the radius measured from the centre 
of the bubble, f the corresponding unit vector and g(r) — 
(r/R) 2 if r < R and g(r) = 1 if r > R. 

The bubble was filled with around 10 4 uniformly dis- 
tributed tracer particles that are advected with the flow. We 
employed a polytropic equation of state with an adiabatic 



index of 7 = 5/3. The simulations with the serial version of 
the code were performed on a SUN ULTRA 10 workstation 
and the parallel version of the code was run on 8 processors 
on an SGI ORIGIN 3000. 

Finally, we should address some issues related to the 
accuracy of these kinds of finite-difference hydrodynami- 
cal simulations. First, one can note that the gas density in 
the bubble increases with time. Any observed diffusion is 
entirely numerical. The boundary between the bubble and 
the ambient medium also becomes less sharp as the simu- 
lation proceeds due to discretisation errors in the advection 
scheme. For a test of the advection algorithm in the ZEUS 
code see Stone & Norman (1992a,b). Simple tests showed 
that during the advection of a sharp discontinuity over a 
grid of 200 zones the discontinuity is spread over 3-4 grid 
cells. Therefore, the smaller features in our simulation are 
more susceptible to advection errors than the larger ones. 

Second, numerical viscosity is responsible for suppress- 
ing small-scale instabilities at the interface between the bub- 
ble and the cooler, surrounding gas. To assess the effects 
of numerical viscosity, we have repeated our simulations in 
2D on grids with 300 2 and 500 2 grid points. Finally, we 
found that numerical resistivity destroys the magnetic field 
upon compression by the shock. With the affordable numer- 
ical resolution the magnetic field orientation reversals are 
pressed in very nearby computational cells. This is discussed 
further in Sec. 3.2. 
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Figure 3. Evolution of the average magnetic field energy density 
~ B 2 , the thermal pressure P, and the electron number density 
n c at the locations of the tracer particles in the sSwB simulation 
(sSwB = strong Shock, weak B-field; see Sect. ^) for notation . 
The impact of the shock compression is visible after around 100 
Myr. Also the rapid decline of the magnetic field strength due to 
numerical resistivity is visible. This decline is not observed in 2-d 
test simulations with either a much higher resolution or a much 
more ordered magnetic field. 



2.2 Radio Maps 

The radio maps are constructed using around 10 4 passively 
advected tracer particles. Initially, each tracer particle is lo- 
cated inside the radio plasma cocoon and is associated with 
the same initial relativistic electron population. Then the 
electron spectrum for each tracer particle is evolved in time 
taking into account synchrotron, inverse Compton, and adi- 
abatic energy losses and gains. In the formalism of Enfilin & 
Gopal-Krishna (2001)^] this spectral evolution requires two 
quantities to be tracked as a function of time. The first is 
the compression history, i.e. 



C(t) = V /V(t), 



(2) 



where Vo and V(t) are the initial and present specific vol- 
umes of the tracer particles. In our calculations we use the 
pressure P(t) at the location of the tracer particles to work 
out the compression via the equation of state: 



C(t) = (P(t)/P ) 1 ^, 



(3) 



where 7 rp is the adiabatic index of radio plasma. Even 
though our MHD simulations have a constant adiabatic in- 
dex of 5/3, we use 7 rp = 4/3 in Eq. ^ in order to simulate the 
higher energy gains of the electrons in the more compressible 
relativistic radio plasma. 

The second quantity that has to be tracked in order to 
evolve the electron spectrum is the energy (or momentum) 
of an electron that was initially injected with infinite energy. 
This cutoff momentum p* of the electron spectrum (here and 

t Equivalent formalisms were derived and applied to numerical 
simulations by other authors, e.g. Matthews & Scheuer ( 1990a, b) 
, and Churazov et al. (2001) . 



Figure 4. Evolution of the radio spectrum of the sSwB sim- 
ulation. The numbers (00, 10, 90) label different simulation 
data-dump time-steps. A single time step is approximately 7 Myr 
long. Clearly visible is the decrease of the cutoff frequency with 
time. Only in the interval 10-20 (70-140 Myr) the shock com- 
pression is able to reverse this. During this phase also the overall 
cmissivity at lower radio frequencies increases mainly due to field 
strength amplification. At all other time-steps the low frequency 
flux decreases artificially, due to magnetic field annihilation by 
numerical field diffusion. 



in the following dimensionless units for electron momenta 
p — P c /(m e c) are used, where m a is the electron rest-mass, 
and c the speed of light) is given by 



^=a Jdt>(u B (t>) + uc(t')) 



(4) 



Here ao = | <tt/ ( m e c) with ctt the Thomson cross-section, 
and ub and uc are the magnetic and cosmic microwave 
background energy densities, respectively. Hence the energy 
evolution of any electron can be calculated using 

Pa 



p(po,t) = 



C(t)-*+Po/P.(*) 



(5) 



The electron spectrum per dimensionless momentum of 
a tracer particle /(p, t) dp for p < p*(t) is given by 



/(p,*) = /o(po(p,*)) 



dpo(p,t) 



(6) 



dp 
where 

po(p,t)=pc(tr*/(i-p/p«(*)). (7) 

If the original distribution function is a power-law 

/o(po) = goPo"" (8) 

for Pmino < Po < PmaxO the resulting spectrum is 

/(p,t)= 5 oC*(t)^V Qc (l-p/p,(i)r~ 2 (9) 

for p m in(t) =p(PminO,i) <P <Pmax(t) =p(p m ax0,t)- 

At fixed intervals during the simulation (up to 100) the 
positions of the tracer particles as well as the gas and mag- 
netic field properties were written out. These were then used 
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to tabulate C(t) and p*(t) interpolating linearly in the inte- 
gration in Eq. m A second programme read in these tables 
to calculate the Stokes I-, Q-, and U-polarisation radio maps 
using the standard integration kernels of synchrotron radi- 
ation theory (Rybicki & Lightman 1979). At this step we 
have specified the parameters of the initial electron spec- 
trum, viewing angle, distance and observing frequency. The 
maps were stored in the FITS format, and visualised with 
the radio-astronomical data reduction package AIPS. 



3 RESULTS 

In the following, results from four simulation runs are pre- 
sented: two simulations with a weak shock (external com- 
pression factor C gas = 2), and two with a stronger shock 
(Cgas = 3.3). For each shock strength, simulations with dy- 
namically unimportant ('weak') and dynamically important 
('strong') field strengths are discussed. In those simulations 
with strong magnetic fields, the magnetic pressure initially 
exceeds the thermal pre-shock pressure. But a short violent 
expansion quickly establishes approximate pressure equilib- 
rium. In the following, we name the simulation with the 
weak shock and the weak field strength 'wSwB', with the 
weak shock and the strong field strength 'wSsB', and like- 
wise. 

Physical units are assumed such that the simulated ra- 
dio relic resembles one which is located in the denser regions 
of a small galaxy cluster. In principle, we could assume phys- 
ical units that correspond to more peripheral relics that re- 
side in lower pressure environments. But i n th is case the arti- 
ficially decaying magnetic fields (see Sec. 3.2) have a greater 
effect on the radio luminosity because the inverse Compton 
losses are bigger than the synchrotron losses. 

The pre-shock external gas density is set to about 5 • 
10 -4 electrons/cm 3 and the temperature to 1-2 keV. The 
simulation box is assumed to have a size of 1 Mpc 3 and to 
be located at a distance of 100 Mpc from the observer. 



GHOST IPOL 100 MHz R = 2 weak B face on view 



3.1 Torus Formation 

The passage of a radio cocoon through a shock wave can be 
seen in Fig. ^| The formation of a torus can be explained 
as follows: At a shock wave, the ram-pressure of the pre- 
shock gas (plus its small thermal pressure) and the thermal 
pressure of the post-shock gas are in balance. But when the 
lighter radio plasma comes into contact with the shock sur- 
face, the ram pressure is reduced at the point of contact. 
The post-shock gas starts to expand into the volume occu- 
pied by the still low pressure radio plasma. This can be seen 
in the second and third displayed panel of Fig. ^| Finally, 
the radio plasma is disrupted and a torus forms. 

A possible evolutionary scenario of the radio morphol- 

fy of the cocoon could be similar to what is shown in Fig. 
In the early stages of the compression the torus has not 
formed and one expects a sheet like radio relic to be visible. 
The relic should get more and more edge brightened while 
the radio plasma accumulates in the torus. Finally, when 
the thin sheet is disrupted, only a torus-like object should 
be visible. 

Thus we expect that the early, sheet-like radio relics 
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Figure 5. Radio emission of the radio bubble in the model wSwB 
before shock passage. The polarisation E- vectors are displayed by 
dashes with the length proportional to the relative polarisation. 
Here and in the following radio maps the flux is given in Jansky 
per simulation pixel (with linear size of 20.6 arcsec) 
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Figure 6. Same as Fig. ^. Late stage of the shock passage in the 
same model wSwB. 



have a 'younger', less steepened spectrum, while the later- 
stage, filamentary or toroidal relics have a much higher spec- 
tral age with a bent, steep radio spectrum. 

These predictions are supported by the fact that the rel- 
atively sheet-like relics in the cluster Abell 2256 have a flat 
spectrum (Rottgering et al. 1994) whereas all the filamen- 
tary relics in the sample observed by Slee et al. (2001) ex- 
hibit steep, bent and therefore 'aged' spectra. We note that 
this statement needs additional observational confirmation. 
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GHOST IPOL 100.000 MHz R = 3.3 weak B face on 
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Figure 7. Facc-on view of the late stage of the shock passage in 
the model wSsB. For details see Fig. ^. 



Figure 9. Face-on view of the late stage of the shock passage in 
the model sSwB. For details see Fig. H. 



GHOST IPOL 100 MHz R = 2 strong B edge on view 




GHOST IPOL 100 MHz R = 3.3 strong B viewing angle 45 degree 
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Figure 8. Edge-on view of the late stage of the shock passage 
in the model wSsB. For details see Fie 



3.2 Global Properties 

In Fig. ^the evolution of the average magnetic field strength, 
the thermal pressure and the electron number density at the 
location of the tracer particles are shown for one simulation. 
The bubble reaches its maximal compression after 75-200 
Myr. The associated increase in the thermal electron den- 
sity is somewhat higher than predicted by simple adiabatic 
compression. This is probably due to gas electrons that have 
diffused numerically into the radio plasma volume. However, 
this does not affect our radio maps. The spatial distribution 
of the relativistic particles is solely determined by the loca- 



Figure 10. 45 degree view of the late stage of the shock passage 
in the model sSsB. For details see Fig. pi 



tions of tracer particles. The compressional state of the rel- 
ativistic electrons is calculated from the fluid pressure and 
not from the density in order to be less sensitive to numerical 
diffusion. 

Moreover, the decline of the magnetic field strength is 
caused by the numerical resistivity (or magnetic diffusivity) 
of the code. The reason for the strong impact of the nu- 
merical diffusion of the magnetic field on our results lies in 
the thin morphology the radio plasma after the shock pas- 
sage. This brings low density and high density regions, and 
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Figure 11. 100 MHz polarisation scans around the radio source 
in the simulation wSsB at an early stage of the shock passage 
with the viewing angle ranging from face-on (viewing angle = 0°) 
over edge-on (viewing angle = 90°) to viewing the back (viewing 
angle = 180°). The different angles <f> = -45, 0, +45, +90 label 
four different routes of the scans around the relic. 

also regions with differently oriented magnetic fields close 
together, producing strong gradients which are sensitive to 
numerical diffusion. Since the radio plasma is compressed 
into filaments that are only a few cells wide, magnetic diffu- 
sion leads to strong artificial annilihation (or reconnection) 
of magnetic fields. 

Two-dimensional resolution studies indicate that this 
artificial decrease of the magnetic field strength should be- 
come small for a 500 3 cell simulation. Since such a simulation 
is at present not feasible for us we will try to focus on robust 
conclusions that are not severely affected by this inaccuracy. 
However, we have made experiments with low resolution and 
large-scale fields in two dimensions. These experiments re- 
vealed that the magnetic field strength increases drastically 
upon compression and, subsequently, remains strong - as is 
expected in the absence of numerical diffusivity. 

3.3 Spectral Aging 

An increase of a factor ~ 5 in the radio luminosity during 
the compression between 70 to 140 Myr can be seen in Fig. 
^. However, the greater increase expected from the analyt- 
ical model of Enfilin & Gopal-Krishna (2001) could not be 
reproduced in our simulations. This is a result of the decay- 
ing magnetic fields due to numerical resistivity, and not a 
failure of the model. In order not to be too much affected by 
the rapidly evolving spectral cutoff we restrict our analysis 
in the following to a low radio frequency, namely 100 MHz. 

3.4 Radio Morphology 

As the shock wave passes the initially spherical radio co- 
coon (Fig. the cocoon is torn into a filamentary structure 
(Figs. and0). In the simulations with weak mag- 

netic fields the final morphology is toroidal for the strong 
shock (simulation sSwB, Fig. m and shows two tori in the 



Figure 12. Same as Fig. gjj, but for a late stage of the shock 
passage. 
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Figure 13. Same as Fig. |l2| but here the angle between the 
projected normal of the shock wave and the E-polarisation vector 
is displayed. Values near zero mean that the projected magnetic 
fields are mostly aligned with the projected shock orientation. 

case of a weak shock wave (wSwB, Fig. ^). Such a double 
torus seems to be observed in the relic in Abell 85 (Slee et 
al. 2001, and Fig. |l|). The degree of polarisation is relatively 
high everywhere. The electric polarisation vectors tend to 
be perpendicular to the radio filaments indicating aligned 
magnetic field structures within them. The reason for this 
is that the compression is mostly perpendicular to the fila- 
ments since this is how they formed. Therefore, the parallel 
field component is preferentially enhanced. 

But in some few spots (see Fig. ^ the polarisation E- 
vectors are aligned with the radio filaments. We presume 
that these spots coincide with regions where the initial mag- 
netic component parallel to the shock was too weak to be 
amplified above the perpendicular components. Such rever- 
sals are indeed observed in polarisation maps of relics but 
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mult i- frequency observations still have to confirm that they 
are intrinsic, and not induced by foreground Faraday rota- 
tion (R. Slee, private communications). 

In those simulations where the magnetic fields are dy- 
namically important nearly always a more complicated mor- 
phology is produced. In the weak shock simulation (wSsB, 
Figs. [7] and ^) the cocoon does not show any toroidal struc- 
ture after shock passage. But filamentary structures can 
still be identified and the E-polarisation vectors seem to be 
approximately perpendicular to the filaments. The edge-on 
projection shows that the radio emitting volume is not re- 
stricted to the shock plane but also shows kinks and exten- 
tions in the flow direction. Similar looking kinks are observed 
in the radio relics in the cluster Abell 3667. These relics are 
also likely to be seen edge-on as the merger shock lies within 
or close to the plane of the sky (Roettiger 1999) . 

In the simulation with a strong shock and strong fields 
(sSsB) a large torus still appears, but a second smaller one is 
dominating the total radio emission (Fig. |l^). The extreme 
brightness of this feature is likely to be an artifact of an 
unlucky realization of the initial random magnetic field con- 
figuration than a typical feature of a real radio relic, since 
a corresponding bright spot already appeared in the initial 
radio cocoon. 



3.5 Radio Polarisation 

Our simulations indicate that for a sufficiently turbulent ini- 
tial magnetic field geometry the total integrated polarisation 
of a radio relic observed face-on should be small. It is of the 
order of the total initial fractional polarisation of the source. 
In the case of a simulated toroidal relic all polarisation ori- 
entations are roughly equally present and therefore cancel 
each other out in the surface integration. But even in a more 
complex relic, as shown Fig. ^| the integrated polarisation 
cancels out in the face-on view. 

In the edge-on view this is different. Now a preferen- 
tial direction exists: E-vectors tend to be aligned with the 
projected shock normal. 

It is interesting to study how reliably the overall polar- 
isation strength and direction can be used as a measure of 
the 3-dimensional shock orientation. 

Analytic predictions for the polarisation strength of 
shock compressed random magnetic fields were already de- 
rived in Laing (1980) and Enfilin et al. (1998). In these 
derivations it was assumed that the compression is linear 
along the flow axis and that the fields follow completely 
passively. The resulting percentage radio polarisation is then 
given by: 



S + l 



sin 2 8 



(10) 



The term C gas denotes the shock compression factor, s the 
spectral index of the electron population, and S the viewing 
angle. If the magnetic fields dominate the total pressure one 
can allow them to expand to pressure equilibrium after shock 
passage. In this case the polarisation 
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(11) 



(Enfilin et al. 1998). This can be seen from the reference 
curves in Fig. |Ti] and 

Also shown in these figures are viewing angle scans of 
the total polarisation of the radio emission in the simula- 
tion wSsB before (Fig. 0) and after (Fig. Q shock pas- 
sage. These figures show that, in addition to the intrinsic 
polarisation pattern of the turbulent fields, the shock pas- 
sage imprints a strong polarisation signature onto the radio 
plasma. 

It is remarkable that the characteristics of the simu- 
lated radio polarisation are similar to the analytic models of 
shocked magnetic fields. If the amplitude of the polarisation 
curves in our simulations are not too strongly affected by ar- 
tifacts, such as the decaying field strength, or the incorrect 
adiabatic index of the radio plasma, we can conclude that 
the shock acceleration and the revived fossil radio plasma 
are not distinguishable by this. Alternatively, the simplistic 
analytic theory of shocked fields seems to be a good guide 
in relating viewing angle and radio polarisation. 

The reason for the similar polarisation pattern in the 
different models is likely that the magnetic fields in all cases 
are nearly completely aligned with the shock plane. In the 
extreme case of completely shock plane aligned fields, but 
with randomly distributed orientations within that plane, 
the polarisation is given by 



> = 



s + l 



s + | 2 



(12) 



is surprisingly similar to that expected in the weak field case 



This is not too far from the other models. 

Whenever the imprinted polarisation dominates over 
the initial intrinsic one, the direction of the E-vector can 
be used relatively reliably to infer the shock normal pro- 
jected onto the plane of the sky (Fig. |l^). The angle with 
respect to this plane can be estimated roughly from the total 
polarisation. 

Thus, in principle, the 3-dimensional orientation (mod- 
ulo a mirror ambiguity) can be derived from the polarisa- 
tion data only. However, these points will need additional 
investigations before it is applicable to real data. Especially 
the role of the initial field geometry has to be investigated. 
Our simulated radio ghosts showed deviations from the ana- 
lytic predictions on the order of the total initial polarisation 
of the source. Such deviations may significantly reduce the 
accuracy of any polarisation based viewing angle estimate. 
However, we think that according to our model there should 
be a correlation between total polarisation and the viewing 
angle due to the field ordering action of the shock wave. See 
Enfilin et al. (1998) for an early attempt to compare these 
quantities for a small smaple of cluster radio relics. 



3.6 Shock Properties 

On the basis of our simulations we speculate that the ob- 
served dimensions of a cluster radio relic with a toroidal 
shape may be used to get a rough estimate of the shock 
strength. This is based on the observation that in the nu- 
merical simulations the radius R of the spherical cocoon and 
the major radius of the torus after the passage of the shock 
are approximately equal. Since the major (R) and minor (r) 
radii of the torus can be read off approximately from a sensi- 
tive high-resolution radio map, the compression of the radio 
plasma by the shock can be estimated. In the idealised case 
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of an initially spherical and finally toroidal radio cocoon, the 
compression factor is given by 

r _ Sphere _ 4 7rfl 3 /3 2R 2 

\4orus 2ir 2 Rr 2 3nr 2 ' 1 j 

Assuming the radio cocoon to be in pressure equilibrium 
with its environment before and after the shock passage, 
the pressure jump in the shock is given by P2/P1 = C 7lp . 
If the adiabatic index of the radio plasma 7rp is assumed to 
be known, e.g 7 rp = 4/3 for an ultra-relativistic equation of 
state, the shock strength can be estimated. But even if this 
assumption is not justified and the geometry deviates from 
the idealised geometries assumed here, the strength of the 
shock wave should be correlated to the ratio of the global 
diameter of a toroidal relic and the thickness of its filaments. 
Unfortunately, the quality of the best current radio maps of 
relics do not yet allow a quantitative comparison of the shock 
strength by comparing the R/r ratios of toroidal relics. But 
these maps demonstrate that the necessary sensitivity and 
resolution might be reached soon. 

If, furthermore, the strength of the shock wave of well 
resolved cluster radio relics can be estimated independently 
from X-ray maps of the IGM, it would be possible to directly 
measure the adiabatic index of radio plasma. This can be 
done by estimating the slope of a log(P2/Pi) versus log(7?/r) 
plot for a sample of well observed radio relics and their shock 
waves: 

log(P 2 /P 1 ) = 2 7 r P log(P/r) - log(37re/2) . (14) 

Here e parametrises the influence of the deviation from the 
above assumed ideal geometry. The values of e should scatter 
from relic to relic, but any systematic correlation of e with 
the shock strength should be weak. Even though the present 
radio and X-ray data do not have the required accuracy yet, 
at some point in the future this method may enable us to 
measure the unknown equation of state of radio plasma. 

In our simulations the adiabatic index of the simulated 
fluid (7r P ,sim = 7 g as = 5/3) and the shock strength are 
known (P2/P1 = 3.5 and P2/P1 ~ 17.4 for the shock com- 
pression factor Cshock = 2 and C s hock = 3.3 respectively). 
Thus we find that the ratio of the length scales R/r w 3 for 
Cshock = 2 and R/r ps 5 for C s h oc k = 3.3. This is roughly 
consistent with our synthetic radio maps. At least the qual- 
itative correlation of shock strength and diameter ratio is 
clearly observed, as can be seen in Figs. |, |, and 0. Even 
in the non-toroidal, filamentary case wSsB the ratio of the 
typical filament diameter to the relic diameter seems to be 
similar to that in the toroidal case, as a comparison of Fig. 
Q with Fig. [] shows. 



4 CONCLUSION 

We have presented 3-D MHD simulations of a hot, magne- 
tised bubble that traverses a shock wave in a much colder 
and denser environment. This is assumed to be a fair model 
for a blob of radio plasma in the IGM which is passed 
by a cluster merger shock wave. We have calculated ra- 
dio polarisation maps for the relativistic electron popula- 
tion and computed the spectrum subject to synchrotron-, 
inverse Compton- and adiabatic energy losses and gains. 



These maps show that the shock wave produces filamen- 
tary radio emitting structures and, in many cases, toroidal 
structures. Such filaments and tori appear to be observed by 
very recent high-resolution radio maps of cluster radio relics 
(Slee et al. 2001). Our simulations find polarisation patterns 
which indicate that the magnetic fields are mostly aligned 
with the direction of the filaments. This also seems to be the 
case for the observed cluster radio relics. 

We argued that the formation of filaments and tori is 
a generic feature of a shock processed hot bubble of plasma 
whose internal sound speed well above the shock speed. We 
expect this to be a robust result that will also be found in 
more realistic simulations, e.g. with the proper equation of 
state, or more realistic initial magnetic field configurations. 

Therefore, we conclude that we have found support for 
the hypothesis that cluster radio relics indeed consist of fossil 
radio plasma that has been compressed adiabatically by a 
shock wave, as proposed by Enfilin & Gopal-Krishna (2001). 
Thus, the historical name 'cluster radio relic' seems to be, 
accidentally or intuitively, an appropriate choice. 

From our simulations several predictions about the 
properties of cluster radio relics can be made - if our sce- 
nario is correct : 

The formation of the tori and filaments is not instan- 
taneous. First, the simulations show a phase in which the 
radio plasma is strongly flattened by the shock. During this 
phase a sheet-like radio relic with a flat spectrum is ob- 
served. Later, the radio plasma moves towards the edges of 
this sheet and finally becomes a torus (or a more compli- 
cated, filamentary structure). Since spectral ageing is likely 
to have affected these later stages, we expect that the fil- 
amentary relics have a steeper, more bent radio spectrum 
than the sheet-like ones. 

Our simulations indicate that the diameter D of the 
bubble of radio plasma remains approximately constant dur- 
ing the passage of the shock. It was also found that the final 
structure consists of radio filaments of small diameter d that 
are distributed (often in form of a torus) in an area of di- 
ameter D. The compression factor of the radio plasma is 
proportional to (D/d) 2 , and this ratio is a measure of the 
shock strength. Thus, the approximate compression factor 
of the radio plasma can in principle be read off a sensitive 
high-resolution radio map. If this number can be estimated 
for a sufficient number of late-stage relics, and if the pres- 
sure jump P2/P1 of the shock wave can be obtained from 
detailed X-ray observations, the unknown adiabatic index 
7rp of the radio plasma can be computed using the relation: 

T^y-d) ■ (15) 

The local radio polarisation reflects the complicated 
magnetic field structures. But for sufficiently turbulent ini- 
tial magnetic fields (as in our simulations) the total inte- 
grated polarisation of a relic reveals the 3-dimensional ori- 
entation of the shock wave. Since the compression aligns the 
fields with the shock plane, the sky-projected field distribu- 
tion is aligned with the intersection of the shock plane and 
the sky plane. Thus, the direction of the total E-polarisation 
vector yields the sky-projected normal of the shock wave. In 
principle the angle between the normal of the shock and 
the plane of the sky can be estimated from the fractional 
polarisation of the integrated flux. Our simulations indicate 
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that the viewing angle dependence is similar to that which 
one would expect if the magnetic fields were embedded in 
the external gas and passively compressed by a shock, even 
if this may not be the physical mechanism. But this state- 
ment should be treated with some caution since the simu- 
lated polarisation may be sensitive to artifacts such as the 
high numerical magnetic diffusivity and the incorrect adi- 
abatic index of the simulated radio plasma (7 rp , s i m = 5/3 
instead of the more likely 7 rp = 4/3). Furthermore, any ini- 
tial global field orientation may survive the compression and 
distort this relation between viewing angle and polarisation 
strength. However, even in this case one expects still a sta- 
tistical correlation between these quantities. These points 
need further investigation. 

We conclude that this work supports the hypothesis 
that cluster radio relics are revived bubbles of fossil radio 
plasma, the so called radio ghosts. Spectral aging arguments 
(Enfilin & Gopal-Krishna 200 f) predict the existence of a siz- 
able population of yet undetected cluster radio relics which 
are only observable with sensitive low frequency radio tele- 
scopes. 
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